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I. INTRODUCTION 



The production of particles with high transverse momenta from the coUision of hadrons and 
nuclei at sufficiently high energies has relevance to constituent interactions at small scales. In this 
regime, it is interesting to search for new physical phenomena in elementary processes such as quark 
compositeness [1], extra dimensions [2], black holes [3], fractal space-time [4] etc. Other aspects 
of high energy interactions arc connected with small momenta of secondary particles and high 
multiplicities. In this regime collective phenomena of particle production take place. The search 
for new physics in both regions is one of the main goals of investigations at Relativistic Heavy 
Ion ColKder (RHIC) at BNL and Large Hadron Collider (LHC) at CERN [5]. Experimental data 
on particle production at high energy and multiplicity provide constraints for different theoretical 
models. Processes with high transverse momenta of produced particles are suitable for a precise test 
of perturbative Quantum Chromodynamics (QCD). The soft regime is preferred for verification of 
non-perturbative QCD and investigation of phase transitions in non-Abelian theories. 

One of the methods allowing systematic analysis of data on inclusive cross sections over a wide 
range of the collision energies, multiplicity densities, transverse momenta, and angles of the pro- 
duced particles is based on the ^;-scaling observed in high energy proton-(anti)proton collisions (see 
Ref. [6] and references therein). The approach to the description of the inclusive spectra reflects 
the principles of locality, self-similarity, and fractality of hadron interactions at a constituent level. 
It takes into account the structure of the colliding objects, interaction of their constituents, and 
processes of particle formation over a wide scale range. The analyzed data include processes which 
take place at nucleon scales as well as at small scales down to 10~^ Fm. The presentation of experi- 
mental data in this approach is given in terms of the scaling function t}){z) and the scaling variable z. 
Both are constructed using the kinematical variables, the experimentally measured inclusive cross 
section Ecfia/dp^, the multiplicity density dN/drj, and some model parameters which allow physical 
interpretation. 

It was shown [6] that the scaling behavior of tp{z) is valid for different types of the produced 
hadrons. The scaling function demonstrates two regimes. The first one is observed in the high- 
and the second one in the low-pT region. The hard part of the inclusive spectra is described by the 
power law, tp{z) ~ z~^, with a constant value of (3 at large z. The self-similar features of particle 
production dictated by the 2;-scaling at large z give strong restriction on the asymptotic behavior of 
the spectra in the high-p^ region. This provides suitable constraints on phenomenology of parton 
distribution functions and fragmentation hmctions which are needed to verify the perturbative QCD 
in hadron collisions. The soft regime of particle production demonstrates flattening of the scaling 
function tp{z) for small z. The behavior of tp{z) at small z is of interest to study the non-perturbative 
QCD in the processes with low-pT- We consider that this region is most preferable in searching for 
phase transitions and study of collective phenomena in multiple particle systems. 

In this paper we show that in the high energy pp- and pp-coUisions the shape of the scaling 
function at low z is independent of the type of the inclusive hadron including production of the 
hadrons with heavy flavor content. A saturation of ip{z) with decreasing z is observed. The single 
parameter c which controls the behavior of ^{z) at low z is interpreted as a specific heat of the 
produced medium. The scaling in pp- and pp-collisions is consistent with a constant value of c. 
Search for a possible change in this parameter is of interest especially for soft processes with high 
multiplicities. Such a change could be an indication of a phase transition in the matter produced 
in high energy collisions of both hadrons and nuclei. 

The paper is organized as follows. A concept of the ^;-scaling and the method of construction 
of the scaling function ^p{z) are briefly described in Sec. II. The properties of energy, angular, and 
multiplicity independence of the scaling function are mentioned in Sec. III. The flavor independence 
and saturation of ip{z) at low z are demonstrated in Sec. IV. A microscopic picture of the constituent 
subprocess is analyzed in Sec. V. A relation of the scaling variable z to some thermodynamical 
quantities is discussed in Sec. VI. Here we consider possible manifestations of phase transitions 
and their effects to the parameters used in the z-presentation of inclusive spectra. Conclusions are 
summarized in Sec. VII. 
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II. Z-SCALING 



In this paper we follow the version of the i;-scaling presented in Ref . [6] . Let us briefly remind 
the basic ideas of this concept. At sufficiently high energies, the collision of extended objects like 
hadrons and nuclei is considered as an ensemble of individual interactions of their constituents. The 
constituents are partons in the parton model or quarks and gluons in the theory of QCD. A single 
interaction of constituents is illustrated in Fig. 1. 




FIG. 1. Diagram of the constituent subprocess. 

Structures of the colliding objects arc characterized by parameters Si and 62- The constituents 
of the incoming hadrons (or nuclei) with masses Mi , M2 and momenta Pi , P2 carry their fractions 
Xi,X2- The inclusive particle carries the momentum fraction ?/„ of the scattered constituent with 
a fragmentation characterized by a parameter Eq. A fragmentation of the recoil constituent is 
described by £(, and the momentum fraction yj,. Multiple interactions are considered to be similar. 
This property represents a self-similarity of the hadronic interactions at the constituent level. 



A. Momentum fractions xi,X2,ya, and yb 

The idea of the ^;-scaling is based on the assumption [7] that gross features of an inclusive particle 
distribution of the reaction 

Mi+M2^mi + X (1) 

can be described at high energies in terms of the kinematical characteristics of the corresponding 
constituent subprocess. We consider the subprocess to be a binary collision 

(a;iMi) + {X2M2) {mi/ya) + {xiMi + X2M2 + m2/yb) (2) 

of the constituents {xiMi) and {X2M2) resulting in the scattered [mi/ya) and recoil (X1M1+X2M2 + 
ff^ilVb) objects in the final state. The produced secondary objects transform into real particles after 
the constituent collisions. The registered particle with the mass mi and the 4-momentum p and its 
hadron counterpart, moving in the opposite direction, carry the momentum fractions j/a and t/t, of 
the scattered and recoil systems, respectively. The momentum conservation law of the constituent 
subprocess is connected with a recoil mass which we write in the form 

{xiPi + X2P2 - p/Vaf = {xiMi + X2M2 + ma/ybf. (3) 

The associate production of (7712) ensures conservation of the additive quantum numbers. Equation 
(3) is an expression of the locality of the hadron interaction at a constituent level. It represents a 
kinematical constraint on the momentum fractions Xi, X2, ya, and yb which determine a subpro- 
cess (2). 
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Structure of the colliding objects and fragmentation of the systems formed in the scattered and 
recoil directions are characterized by the parameters 61,62 and Soi Cfci respectively. We connect the 
structural parameters with the corresponding momentum fractions by the function 

nixi,X2,ya, Vb) = (1 - X,Y' (1 - X2f' (1 - VaY-il ~ VbY' ■ (4) 

Physical interpretation of Q. is given by its proportionality to relative number of all such constituent 
configurations in the reaction (1) which contain the configuration defined by the fractions Xi,X2,ya, 
and Ub- The function plays the role of a relative volume which occupy these configurations in the 
space of the momentum fractions. It was found that the structural parameters 6i, 62, ea, and Cb have 
constant values at high energies. They are interpreted as fractal dimensions in the corresponding 
space of the momentum fractions. For proton-proton collisions we set 61 = 62 = 6. In the case 
of nucleus- nucleus collisions there are relations 61 — A16 and 62 = A26, whore Ai, A2 are atomic 
numbers [8]. We assume that the fragmentation of the objects moving in the scattered and recoil 
directions can be described by the same parameter Ca = Cb = cp which depends on the type (F) 
of the inclusive particle. For given values of 6 and ep, we determine the fractions xi, X2, ya, and 
yb in a way to maximize the function X2, j/a, J/ft), simultaneously fulfilling condition (3). The 

momentum fractions Xi and X2 obtained in this way can be decomposed as follows 

a;i = Ai+xi, X2 = X2+X2, (5) 

where Ai_2 = Ai 2(2/0. yt) and xi.2 = Xi, 2(2/0, 2/5) arc simple specific functions [6] of ya and yb- Using 
the decomposition, the subprocess (2) can be rewritten into a symbolic form 

Xi+X2^ (Ai + A2) + (xi + X2)- (6) 

This relation moans that the A-parts of the interacting constituents contribute to the production of 
the inclusive particle, while the x-parts are responsible for the creation of its recoil. The maximum of 
the function (4) with the condition (3) can be obtained by searching for the unconstrained maximum 
of the function 

F{ya,yb) = ^{xi{ya,yb),x2{ya,yb),ya,yb) (7) 

of two independent variables 2/0 and yb- There exists a single maximum of F{ya,yb) for every 
momentum p in the allowable kinematical region. Values of ya and 2/6 corresponding to this maximum 
have been determined numerically. Having obtained the fractions xi, X2, ya, and yb, we evaluate 
the function O according to (4). For fixed numbers 6 and we obtain in this way the maximal 
value of ri for every momentum p of the inclusive particle. 

Since the momentum fractions are determined by means of the maximization of the expression 
(4), they implicitly depend on 6 and ep- The parameter ep enables to take effectively into account 
also prompt resonances out of which the inclusive particle of a given type may be created. At 
fixed mass parameter m2, larger values of ei? correspond to smaller and yb, which in turn give 
larger ratios m2/yb and mi/ya- In our phenomenological approach this means that production 
of the inclusive particle (rni) and its counterpart (m2) is a result of fragmentation from larger 
masses which mimic in a sense processes with prompt resonances. Values of these parameters are 
determined in accordance with the experiment and are discussed in the next sections. 



B. Scaling variable z and scaling function ip{z) 

The self-similarity of hadron interactions reflects a property that hadron constituents and their 

interactions arc similar. This is connected with dropping of certain dimensional quantities out of 
the description of physical phenomena. The self-similar solutions are constructed in terms of the 
self-similarity parameters. We search for a solution 

, / ^ 1 c^f 
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depending on a single self-similarity variable z. Here oinei is an inelastic cross section of the reaction 
(1) and N is an average particle multiplicity. The variable z depends on momenta and masses of 
the colliding and inclusive particles, structural parameters of the interacting objects, and dynamical 
characteristics of the produced system. We define the variable z as follows 

z = zofl~'^, (9) 

where 



° {dN,h/dv\orm ^ ' 

and fl given by (4). For a given reaction (1), the variable z is proportional to the transverse kinetic 

energy y^sX of the constituent subprocess (2) consumed on the production of the inclusive particle 
(mi) and its counterpart (m2). The energy ,/s± is determined by the formula 

V^ = Ta+n, (11) 

where 

Ta = VaiV^ - MiAi - M2X2) - mi, n = Vbi^ - Mixi - M2X2) - ma. (12) 
The terms 

Vil = [(AlPi + X2P2f]'/', = [(XlPl + X2P2f]'/^ (13) 

represent the energy for production of the secondary objects moving in the scattered and recoil 
direction, respectively. The quantity dNch/dr]\o is the corresponding multiplicity density of charged 
particles in the central region of the reaction (1) at pscudorapidity 77 = 0. The multiplicity density 
in the central interaction region is related to a state of the produced medium. The parameter 
c characterizes properties of this medium [6]. It is determined from multiplicity dependence of 
inclusive spectra. The mass constant m is arbitrary and we fix it at the value of nuclcon mass. 

The scaling function ip{z) is expressed in terms of the experimentally measured inclusive cross 
section Ed'^a/dp^, the multiplicity density dN/dr] at pseudorapidity r], and (Tinei- Exploiting the 
definition (8) one can obtain the expression [6] 

where s is the square of the center-of-mass energy and J is the corresponding Jacobian. The 
multiplicity density dN/dr] in the expression (14) concerns particular hadrons species. It depends 
on the center-of-mass energy, on various multiplicity selection criteria, and also on the production 
angles at which the inclusive spectra were measured. The procedure of obtaining the corresponding 
values of dN/dr} from the pr spectra is described in Ref. [6]. The function ^{z) is normalized as 
follows 

1*00 

ij{z)dz = 1. (15) 







The above relation allows us to interpret the function ip{z) as a probability density to produce an 
inclusive particle with the corresponding value of the variable z. 



III. PROPERTIES OF THE SCALING FUNCTION 



Let us remind the properties of the ^;-presentation of experimental data already found in proton- 
(anti)proton collisions at high energies. These are the energy, angular, and multiplicity independence 
of the scaling function tp{z) for different types of hadrons, direct photons, and jets confirmed by 
numerous data obtained at ISR, SppS, Tevatron, and RHIC. 
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A. Energy independence of tp{z) 

The energy independence of the ^-presentation of inclusive spectra means that the shape of 
the scahng function is independent on the coUision energy y/s over a wide range of the transverse 
momentum pt of the inclusive particle. Results on the energy independence of the 2;-scaling for 
hadron production in proton-proton collisions were presented in Ref. [6]. The analyzed data [9 17] 
include negative pions, kaons, and antiprotons measured at FNAL, ISR, and RHIC energies. The 
spectra were measured over a wide transverse momentum range j3t = 0.1 — 10 GeV/c. The cross 
sections decrease from 10^ to 10"^'^ mb/GeV^ in this range. The strong dependence of the spectra on 
the collision energy increases with transverse momentum. The independence of the scaling function 
tp{z) on -y/s was found for the constant values of the parameters c ~ 0.25 and S = 0.5 for all types 
of the analyzed hadrons {■K,K,p). The value oi ep increases with the the mass of the produced 
hadron. None of these parameters depends on the kinematical variables. This was demonstrated for 
charged particles, negative pions, kaons, and antiprotons in the range of -/s = 19 — 200 GeV. The z- 
presentation of particle spectra produced in proton-antiproton collisions was studied as well. In this 
case the energy independence of the scaling function was obtained for different set of parameters. 
The scaling functions for proton-(anti)proton collisions have similar shapes but differ in the region 
of large z. The different values of the slope parameter (3 (f)pp > /3pp) of the scaling function, 
tp{z) ~ , for charged hadron, direct photon, 7r°-meson, and jet production in pp and pp collisions 
at high z were found in Ref. [18]. 

B. Angular independence of ^^{z) 

The angular independence of the z-presentation of inclusive spectra means that the shape of the 
scaling function is independent of the angle 9cms of the produced particles over a wide range of 
the transverse momentum px- Results on the angular properties of the ^;-scaling in proton-proton 
collisions are presented in Ref. [6]. The analyzed experimental data [19] on angular dependence 
of inclusive spectra include negative pions, kaons, and antiprotons measured at ISR energies. The 
angles cover the range 6 cms = 3° — 90°. The central and fragmentation regions are distinguished by 
a different behavior of differential cross sections. The analysis included the transverse momentum 
spectra of charged hadrons for 9cms = S'' and 13" [20] obtained by the BRAHMS Collaboration 
at RHIC. The data cover a wide range of the transverse momenta pT = 0.25 — 3.45 GeV/c of 
the produced hadrons at the collision energy y/s = 200 GeV. The angular independence of the 
scaling function ■tp{z) for charged particles, negative pions, kaons, and antiprotons was obtained 
for the same values of the parameters c, 5, and ep which give the energy independence of ip{z). 
This represents different values of cp for pions, kaons,... and the same values of 5 and c for all of 
them. The scaling function is sensitive to the value of m2 for small production angles 6cms- This 
parameter was determined from the corresponding exclusive reactions at the kinematical limit (for 
Xi = X2 = Va = Vb = !)• Using Eq. (3), this gives m2 = m{'K~^), = m{K^), and m2 = m{p), for 
the inclusive production of tt", K~ , and antiprotons, respectively. Note that the charged hadron 
multiplicity density dNchld'q\r)=a represents an angular independent factor in the definition of the 
variable z which is the same for all particle species. On the contrary, the scaling function for pions, 
kaons, antiprotons,... is normalized (14) to the angular dependent multiplicity density dN/dr] of the 
corresponding particles, respectively. 

C. Multiplicity independence of V'(^) 

The multiplicity independence of the ^-presentation of inclusive spectra means that the shape of 
the scaling function tp{z) does not depend on the multiplicity selection criteria characterized by the 
different values of dNch/drj. The multiplicity density influences the shape of the inclusive spectra 
especially at high transverse momenta [21]. At low px the dependence of the cross sections on 
dNch/drj concerns mainly the absolute values and much less the shape. Results on the multiplicity 
independence of the ^;-scaling in proton-proton collisions were presented in Ref. [6] . The analysis was 
performed using the /f^-meson and A-baryon spectra [22] obtained by the STAR Collaboration at 
\fs = 200 GeV for different multiplicity classes. The charged multiplicity density was varied in the 
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range dNch/dv = 1.3 — 9.0. The transverse momentum distributions were measured in the central 
rapidity range \ri\ < 0.5 up to the momentum px = 4.5 GeV/c. The multiplicity independence of 
the scaling function provides a strong restriction on the parameter c. Both data favor the same 
vahie of c = 0.25 as was obtained from the energy independence of the z-scaling in proton-proton 
collisions. Similar applies to the px distributions of charged particles [21] associated with the 
different multiplicity criteria. 

An analysis of the spectra of direct photons and jets gives somewhat different results in the sense, 
that the energy, angular, and multiplicity scaling was obtained by other set of the parameters c, 6, 
and ep- Nevertheless, the corresponding parameters do not depend on the kinematical variables, 
similarly as in the case of the ^-scaling for identified hadrons produced in proton-(anti)proton 
collisions. 



IV. NEW PROPERTIES OF Z-SCALING 



In this section we study possibility of a unified description of the particle spectra of different 
hadrons using properties of their ^-presentation. The analysis is based on the observation that 
simultaneous energy, angular, and multiplicity independence of the ^-scaling for negative pions, 
kaons, and antiprotons produced in proton-proton collisions gives the same shape of the scaling 
function ^jj{z). This flavor independence of tp{z) is confirmed here for other inclusive particles 
including the particles with heavy quarks. The independence in proton-antiproton interactions is 
observed with the exception of large z, where the functions ipi^) for PP~ ^-nd pp-collisions mutually 
differ. The transverse momentum spectra of J/tp [23] and T [24] mesons measured at the Tevatron 
energies ^/s = 1800 and 1960 GeV make it possible to investigate the behavior of the scaling function 
in the region of very small z (up to 10~^). In this region we observe a saturation of ip{z) which can 
be approximated by a constant. The saturation and flavor independence of the scaling function for 
different hadrons with light and heavy quarks are confirmation of the factorization 

i^i^ da 

i,{z)— (16) 



dzdr) drj 

of the inclusive cross sections. This property of the z-scaling is valid for pp-collisions in a large 
range of kinematical variables. The factorization of the differential cross section in the variable z 
and the pseudorapidity 77 was demonstrated [6] as the angular independence of the 2;-scaling for 
Ocms = 3'^ — 90*^. The scaling property at small angles gives strong restriction on the parameter m2 
which was found to be equal to the mass of the inclusive particle for the negative pions, kaons, and 
antiprotons. The same relation, m2 = mi , is used here for all types of the inclusive hadrons. 



A. Flavor independence of V'(^) 

Evidence of flavor independence of z-scaling was for the first time noted in Ref. [18]. It was 
found that the value of the slope parameter /3 of the scaling function ip{z) at high z is the same for 
different types of produced hadrons (tt, K,p). The hypothesis was later supported by the results of 
an analysis of hadron {-k"^'^ ,K,p) spectra for high pT in pp- and p^-coUisions. Here we show that 
flavor independence of the z-presentation of hadron spectra is valid for different hadrons over a wide 
range of the variable z. Hence we exploit the scaling transformation 

z — > apz, tp ~* cep^tp (17) 

for comparison of the shape of the scaling fimction ijj{z) for different hadron species. The parameter 
ap is a, scale independent quantity. The transformation does not change the shape of ipi^)- It 
preserves the normalization equation (15) and does not destroy the energy, angular, and multiplicity 
independence of the z-presentation of particle spectra. 

Figure 2(a) shows the z-presentation of the spectra of negative pions, kaons, antiprotons, and 
A's produced in pp-collisions over the range y/s = 19 — 200 GeV and 9cms = 3° — 90°. The symbols 
represent data on differential cross sections measured in the central [9-11,13,15,16,25] and fragmen- 
tation [19] regions, respectively. The analysis comprises the inclusive spectra of particles [26,27] 
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measured up to very small transverse momenta {pr — 45 MeV/c for pions and pt — 120 MeV/c 
for kaons or antiprotons) . One can see that the distributions of different hadrons are sufficiently 
well described by a single curve over a wide z-range (0.01 — 30). The function 'il'{z) changes more 
than twelve orders of magnitude. The solid lines represent the same curve shifted by multiplicative 
factors for reasons of clarity. The same holds for the corresponding data shown with the different 
symbols. The indicated values of the parameter ep (ctt = 0.2, ex — 0.3, ep ~ 0.35, ea — 0.4) are con- 
sistent with the energy, angular, and multiplicity independence of the z-presentation of spectra for 
different hadrons. The parameters were found to be independent of kinematical variables {^/s,pt, 
and 9 cms)- The scale factors ap are constants which allow us to describe the z-presentation for 
different hadron species by a single curve. The estimated errors of are at the level of 20%. 




(a) (b) 

FIG. 2. The flavor independence of z-scaling. The spectra of 7r~, K~ , p, A (a) and p, uj, 
4>, K* , S (b) hadrons produced in pp-collisions in z-presentation. Data are taken from Refs. 
[9-11,13,15,16,19,25-32]. The solid line is obtained by fitting the data. 



Figure 2(b) shows similar results for other hadrons (p, cj, 0, -ftT*, S) produced in pp-coUisions at 
\/s = 200 GeV and 9cms = 90°. The experimental data [28-32] on inclusive spectra are compared 
with the pion distributions [15] measured at RHIC. The shape oftp{z) for these particles is described 
by the same curve (solid line) as depicted in Fig. 2(a). The black circle at the lowest z ~ 0.007 
corresponds to the STAR data [31] on K* resonances measured in the region where the scaling 
function is saturated. Based on the obtained results we conclude that RHIC data on pp-collisions 
confirm the flavor independence of the z-scaling including the production of particles with very 
small Pt- 

The inclusive spectra of heavier hadrons {J/ip, , B, T) [23,24,33] obtained at the Tevatron en- 
ergies y/s =1800 and 1960 GeV allow us to verify the new property of the z-scaling in pp-collisions. 
The data include measurements up to small transverse momenta (pT — 125 MeV/c for charmonia, 
Pt — 290 McV/c for bottomia, and pt — 500 MeV/c for S-mesons). Figure 3(a) shows the trans- 
verse momentum spectra of J/iPtD^,B, and T mesons in the z-presentation. The scaling function 
is the same for hadrons with light and heavy flavors produced in pp- and pj>-collisions in the range 
z = 0.001 — 4. This is indicated by the same line as in Fig. 2. The corresponding values of the 
parameters ap and ep are found to be aj/^ = 0.23, ap)0 ~ 0.23, as — 0.12, q;x(is) — 0.15 and 
ej/i/; = CD" — 0.4, en — 0.4, eT(is) — 0.4, respectively. Figure 3(b) demonstrates results of 
combined analysis of the RHIC [34] and Tevatron [23] data on J/ip-meson spectra measured in pp- 
and pp-collisions at different energies -^/s = 200, 1800, 1960 GeV and angles 6cms = 22°, 90" in the 
z-presentation. The solid line is the same one as shown in Figs. 2 and 3(a). 
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FIG. 3. The flavor independence of z-scaling. Tfie spectra of J/tp,D'^,B,T, tt^ (a) and J/tp (b) 
mesons produced in pp / pp-co\\isions in z-presentation. Data are taken from Rcfs. [23,24,33,34]. Tfie 
solid line is the same as shown in Fig. 2. 



From the performed analysis we conclude that ISR, RHIC, and Tevatron data on inclusive spectra 
manifest the flavor independence of the scaling function 'ijj{z) over a wide range of z. We would like 
to stress that the obtained result is based on pT distributions of the cross sections Ecfia/dp^ which 
reveal strong dependence on the energy, angle, multiplicity, and type of the produced particle. 



B. Saturation of ^{z) 

The description of particle spectra in terms of the variable z depends on the parameters Si^2, £f, 
and c. In the region of high z, the production of different hadrons is characterized by different values 
oi ep- A reliable estimation of this parameter requires processes with high pT- The same applies 
to Si^2 which relates to the structure of the colliding particles in the initial state. On the contrary, 
the low-z region corresponds to the small values of pt and low momentum fractions xi, X2, Ua, and 
yt- In this region z tends to zero as pr decreases or ^/s increases. Here is f2 ~ 1 and the scaling 

1 /2 

variable can be approximated by 2 ~ Sj^ / {dN cu / d'n\oY ■ The behavior of il}{z) at low z is therefore 
governed by the single parameter c. 

As seen from Figs. 2 and 3, the z-presentation of hadron distributions demonstrates weak 
dependence on the variable z in the soft (low px) region. The data on the pion, kaon, antiproton 
and especially on J and T spectra manifest saturation in the range z = — 10~^. This regime 
is characterized by an approximate constant behavior of the scaling function, 7/^(2) ~ const. The 
similar dependence is observed for other hadrons {K* , B, p,..) at low pT- A characteristic of the 
saturation is the slope P of the scaling function which is diminishing with the decreasing z. The 
value of (3 is approximately zero for z — 10^^ — 10^^. 

One can assume that the asymptotic behavior ofip{z) at z — > is universal and reflects properties 
of the produced system consisting of its constituents (hadrons or quark and gluons). The universal 
scaling behavior in this region suggests that mechanism of particle production at low px is governed 
by soft self-similar processes which reveal some kind of a mutual equilibrium leading to the observed 
saturation. 
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V. A MICROSCOPIC PICTURE OF CONSTITUENT SUBPROCESSES 



The approach based on the ^-scaling concept allows us to develop a microscopic scenario of 
particle production in terms of the constituent interactions. Here we discuss some features of this 
scenario. We attribute the quantity z to any inclusive particle in the reaction (1). The scaUng 
variable z has character of a fractal measure. It consists of the finite part z^ and of the divergent 
factor The factor describes a resolution at which an imdcrlying subproccss can be singled 
out of the inclusive reaction (1). The ^{x\,x-2, Ha, Vb) is proportional to number of all configurations 
containing the incoming constituents which carry the fractions Xi and X2 of the momenta Pi and P2 
and which fragment to the inclusive particle (mi) and its counterpart (7712) with the corresponding 
momentum fractions ya and yi,. The parameters 61 and 62 have relation to fractal structure of the 
colliding objects (hadrons or nuclei). They are interpreted as fractal dimensions of the colliding 
objects in the space of the momentum fractions. The parameters and €1, characterize the fractal 
behavior of the fragmentation process in the final state. A common property of fractal measures 
is their divergence with the increasing resolution. The scaling variable has this property as the 
resolution f]~^ goes to infinity, 

z{n) ^ 00 , if fi-i ^ 00. (18) 

It means that z is the scale dependent quantity. For an infinite resolution the whole reaction (1) 
degenerates to a single subprocess (2), all momentum fractions become unity (a;i = X2 = Va = 
yb = 1) and = 0. This kinematical limit corresponds to the fractal limit 2; = 00. In the general 
case, the momentum fractions are determined from a principle of a minimal resolution of the fractal 
measure z. The principle states that the resolution should be minimal with respect to all 
binary subprocesses (2) in which the inclusive particle mi with the momentum p can be produced. 
It fixes the values of the corresponding momentum fractions xi,X2,ya, and yi, and singles out the 
most effective binary subprocess which underlies the inclusive reaction (1). 

A. Momentum fractions versus pr and \/s 

The method of determination of the momentum fractions makes it possible to analyze kinematics 
of the constituent interactions in the framework of the developed approach. The study of xi —X2 and 
ya — yt correlations and their dependencies on the collision energy and transverse momentum of the 
inclusive particle gives us possibility to look at microscopic picture of the underlaying subprocesses. 

The fractions Xi and X2 characterize amount of the energy (momentum) of the interacting 
protons (antiprotons) carried by their constituents which undergo the binary collision (2) that 
underlies the inclusive reaction (1). Figure 4(a) shows the dependence of the fraction xi on the 
transverse momentum pt of the negative pions, kaons, and antiprotons produced in pp-coUisions 
at y/s = 19, 53, 200 GeV and 9cms = 90". The fractions xi and X2 are equal each other in that 
case. They increase nearly linearly with the transverse momentum px- For fixed px, the fraction xi 
decreases as the collision energy ^/s increases. The Xi is larger for the production of heavy particles 
as compared with light ones. The kinematical limit of the reaction (1) corresponds to xi = X2 = 1 
at any collision energy and for any type of the inclusive particle. This can be seen in Fig. 4(a) for 
\/s = 19 GeV where the fraction xi approximates unity at pt — 9 GeV/c for all three particles. 

The correlation xi — X2 for the ,K~, and p production at the energy -/s = 53 GeV and for 
various detection angles is shown in Fig. 4(b). The central {Ocms = 90'^, 58°, 40°) and fragmentation 
{Ocms = 2.86°) regions are distinguished by a different mutual behavior of xi and X2. The xi — X2 
correlation at Ocms = 90° is depicted as a strait line. For other angles belonging to the central in- 
teraction region, the correlation becomes more complicated and depends on the particle type. Both 
tendencies diminish with the increasing pT where the fraction X2 demonstrates linear dependence 
on xi over a wide range. This is a common feature for all hadron species. The situation is different 
in the fragmentation region. The fraction X2 is much less than Xi and strongly depends on the type 
of the inclusive particle at the small angle Ocms = 2.86°. The heavier particle the larger X2- The 
increase of X2 with xi is very slow in the range of small angles. This is changed dramatically near 
the kinematical limit {xi ~ 1) where the fraction X2 begins to grow larger. 



10 




0.0 I I 0.0 I I 

2 4 6 8 10 0.0 0.2 0.4 0.6 0.8 1.0 

Pt (GeV/c) X, 

(a) (b) 

FIG. 4. (a) The dependence of the fraction xi on the transverse momentum px at 6cms — 90°. (b) 

The correlation between the fractions xi and X2 at ~ 53 GeV in the central and fragmentation 
regions. 



The dependence of the momentum fractions j/q and yt on the kinematical variables (p-r, ^cms, ^/s) 
describes features of the fragmentation process. The fraction ya characterizes dissipation of the 
energy and momentum of the object produced by the underlying constituent interaction into the 
near side of the inclusive particle. This effectively includes energy losses of the scattered secondary 
partons moving in the direction of the registered particle as well as feed down processes from prompt 
resonances out of which the inclusive particle may be created. The fraction yb governs the recoil mass 
in the constituent subprocess. Its value characterizes the dissipation of the energy and momentum 
in the away side direction of the inclusive particle. 
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FIG. 5. The dependence of the fractions ya (a) and yh (b) on the transverse momentum pT for 
Tr^,K^, and p produced in the pp-coUisions at = 19,53, and 200 GeV in the central rapidity 
region. 
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Figure 5(a) shows the dependence of ya on the transverse momentum pT for the negative pions, 
kaons, and antiprotons produced in pp-colhsions at the energy ^/s = 19, 53, 200 GeV and Ocms = 90". 
An curves demonstrate a non-hnear monotonic growth with pT- It means that the energy dissipation 
associated with the production of a high pT particle is smaher than for the inclusive processes 
with lower transverse momenta. This feature is similar for all inclusive reactions at all energies. 
Decrease of the fractions ya with the increasing collision energy is another property of the considered 
mechanism. It corresponds to more energy dissipation at higher energies. This can be due to the 
larger energy losses and/or due to the heavy prompt resonances. The third characteristic is a slight 
decrease of ya with the mass of the inclusive particle. It implies more energy dissipation for creation 
of heavier hadrons as compared with hadrons with smaller masses. The asymptotic value oi ya — ^ 
is reached at the kinematical limit for all particle species. 

The dependence of t/b on pT reflects kinematical properties of the recoil system. One can see 
from Fig. 5(b) that yb is nearly independent olpT- It is smaller than ya and decreases with the 
increasing collision energy -^i. The values of yb are larger for particles with higher masses. The 
qualitative properties of the px dependence of yb are similar for different hadrons. The small values 
of ya mean that the momentum balance in the production of an inclusive particle from a subprocess 
is more likely compensated with many particles with smaller momenta than by a single particle with 
a higher momentum moving in the opposite direction. 
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FIG. 6. The dependence of the fractions ya and yb on the transverse momentum px for (a) tt~ , 
K~ , K*, p, uj, (j), p. A, S~ hadrons produced in pp-collisions at -/s = 200 GeV and (b) J/'ip, D^, 
B, T(IS') mesons produced in pp-collisions at -y/s = 1800, 1960 GeV. 



The pr-dependence of ya and yb for other hadrons produced in pp-collisions have similar behav- 
ior. This is demonstrated in Fig. 6(a) on the data [15,16,25,28-32] at = 200 GeV and Ocms = 90° 
obtained at RHIC. The hollow symbols demonstrate growth of ya with px- The full symbols show 
flattening of yb in the range of pt = 0.2 — 10 GeV/c. Note that both fractions become equal each 
other, ya — yb, for heavier particles and low transverse momenta. This together with mi — m2 
means that, at low px, the objects produced in the constituent collision into the near- and away- 
side direction have equal masses. The feature is even more pronounced for particles with heavy 
quarks. It is well seen from Fig. 6(b) where the momentum dependence of ya and yb is demon- 
strated for J/tp,D'^,B, and T-mesons produced in pp-collisions at the Tevaron energies y/s = 1800 
and 1960 GeV. The spectra correspond to the central rapidity range. The pr-dependence of ya and 
yb has similar behavior as for data obtained at the RHIC energy. 
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There is, however, a particularity in the absolute values of the momentum fractions which concern 
the J/i/j-meson production. The corresponding curves for ija and yh lay much below than what one 
could expect when compared with other particles, even at high pT- This is a consequence of the 
relatively large value of ej/^ = 1 which follows from the requirement of the energy independence 
of ^'(•z) for the charmonium production (see Fig. 3). In our approach it means that J/'(/i-meson is 
produced from larger masses accompanied with extra large dissipation of energy in the final state. 
This exceptional property predestinate the J/^/;-meson to be a suitable probe in AA-collisions where 
the energy losses can be sensitive to different phases of the created matter. 
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FIG. 7. The dependence of the fractions ya and j/f, on the transverse momentum for J /ip (a) and T 
(b) mesons produced in pp- and pp-collisions at = 200, 1960 GeV and 39, 1800 GeV, respectively. 



A comparison of the momentum fractions ya and yb for the charmonium production in pp [34] 
and pp [23] collisions at ^/s — 200 and 1960 GeV is shown in Fig. 7(a). The fractions follow the 
general trend i.e. a decrease with the collisions energy, growth of ya, and approximate constancy of 
j/f, with pt- The equality of both fractions at low pj- remains preserved independent on y^. Another 
feature typical for all types of the inclusive particles is the angular dependence of the fractions. The 
ya increases significantly with the decreasing angle 9 cms- It means that the energy dissipation from 
the constituent subprocess is much smaller in the fragmentation region as compared to the central 
interaction region. The change of the fraction y\j with the angle 6cms is small. A levelling of j/f, with 
Pt for small angles is observed. 

The sensitivity of the transverse momentum dependence of ya and yb to the mass of the T states 
(15, 25, 35) is demonstrated in Fig. 7(b). The symbols correspond to the experimental data [24,35] 
on cross sections measured at FNAL. It turns out that the momentum fractions change visibly 
with the respective state of the bottomium at smaller collision energy. The heavier mass the larger 
fractions. This is connected with different energy losses which seem to depend on the mass state 
of the T-meson at this energy. At higher energies the momentum fractions become insensitive to 
the single states of this particle. The curves for the 15, 25, and 35 states coincide each other at 
^/s = 1800 GeV in the momentum range pr = 0.5 — 17 GeV/c. 
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B. Recoil mass Mx versus pr and ^/s 



Another characteristic of the constituent interactions is the recoil mass 



Mx ^ xiMi+X2M2 + m2/yb (19) 

released in the underlying subprocess (2). It is defined by the right-hand side of Eq. (3). The 
quantity is proportional to the momentum fractions xi and X2 of the interacting objects with the 
masses Mi and M2. Its relation to the fractions ya and yi, is given by the simple dependencies 
2^1,2 = xi^2iya,yb) listed in Rcf. [6]. 
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FIG. 8. The recoil mass Mx in the constituent subprocess underlying the hadron production in 
pp-collisions at ^/s — 19, 53, and 200 GeV in dependence on the transverse momentum px- 



The method of determination of ya and yb makes them implicitly dependent on S and ep- These 
parameters are obtained from the z-scaling approach to analysis of experimental data on inclusive 
particle distributions. In such a way the recoil mass Mx has an internal connection to the structure 
of the colliding objects, constituent interactions, and process of formation of the individual hadrons. 
In the very vicinity of the kinematical limit, the recoil mass is kinematically bounded to approximate 
the value Mi + AI2 + m2- 

Figure 8(a) shows the dependence of the recoil mass on the transverse momenta of the negative 
pions, kaons, and antiprotons produced in pp-coUisions at the energy ^/s = 19, 53, and 200 GeV in 
the central rapidity region. For the sake of clarity, the values of Mx are presented on a log-scale 
with the multiplication factors 10 and 100 for K~ and p, respectively. All curves demonstrate small 
growth at low pr followed by a successive flattening. They reveal a characteristic increase with 
the collision energy and mass of the inclusive particle. The similar dependencies of Mx on pT for 
other hadrons measured at the RHIC energy ^/s =200 GeV are presented in Fig. 8(b). The data 
correspond to the central interaction region. The dependence of Mx on pT for J/4>,D'^,B, and 
T-meson production in pp-collisions at the Tevatron energies ^/s = 1800 and 1960 GeV is presented 
in Fig. 9(a). 

We would like to draw attention to the large recoil mass for the J/?/'-meson production. It is 
approximately equal to the values of Mx for the T-mesons with a small exception at low transverse 
momenta {px below 7 GeV/c). The mass increases with pT and reaches the value of 95 GeV at 
Pt = 18 GeV/c. The curves for B- and D^-mesons lay considerably below. The extra large values 
of the recoil mass in the charmonium production is connected with the abnormal small values of the 
corresponding momentum fractions ya and yb (see discussion to Fig. 6(b)). A comparison of the 
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momentum dependence of Mx for J/i/^-meson at -^s = 200 and 1960 GeV is shown in Fig. 9(b). 
Here one can see the strong sensitivity of the recoil mass to the colhsion energy ^/s. Its angular 
dependence at ^/s = 200 GeV is rather weak. The approximate constancy of Mx at the angle 
9 cms — 22" corresponds to the levelling of the fraction yi, with pT shown in Fig. 7(a). 
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FIG. 9. The recoil mass Mx in the constituent subprocess underlying production of (a) J/ip, -D", 
B, T(15) mesons in pp-collisions at \/s — 1800, 1960 GeV and (b) J/tJ; mesons in pp/pp- collisions 
at ^/s — 200, 1960 GeV in dependence on the transverse momentum px- 



VI. DISCUSSION 



An analysis of ISR, RHIC, and Tevatron experimental data on inclusive cross sections for different 
hadron species performed in z-presentation shows the flavor independence of the shape of the scaling 
function ip{z). This was demonstrated over a wide range of the kinematical variables for proton- 
proton and proton-antiproton collisions. A saturation of ip{z) was found in the region z < 10^^ 
for hadrons with light and heavy quarks. The data on J/V'- and T- meson production confirm 
the saturation up to very small values of z ~ 10"'^. The observed universality is a manifestation 
of the self-similarity of particle production at a constituent level. The dynamical behavior of the 
interacting system revealed in the distributions of the inclusive particles is described by the function 
ip{z) which demonstrates flattening in the experimentally achievable low z-region (Figs. 2 and 3). 
The slope P{z) of the scaling function diminishes with the decreasing z and becomes zero when the 
scaling function flattens out. It was found that for both asymptotic regimes, at low- and high-z, 
the interacting system is characterized by the equation 

,/ ^ ' = const. (20) 
am z 

The obtained results give strong support that flavor independence of V'(^) for variety of light and 
heavy hadrons and its behavior (20) at both asymptotic limits reflect general trend of hadron 
production in the high energy proton-proton and proton-antiproton collisions over a wide range of 
transverse momenta. 
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A. Variable z and entropy 



In this part we discuss analogy between the parameters used in the ^-scaling concept and some 
thermodynamical quantities (entropy, specific heat, chemical potential) of a multiple particle system. 
There exists a connection between the variable z and entropy. The scaling variable is proportional 
to the ratio 



of the transverse kinetic energy ^JsZ and the maximal value of the function 

W{xx,X2,ya,yh) = {dNchldri\Qy Cl{xi,X2,ya,yb) (22) 

in the space of the momentum fractions. The quantity W is proportional to all parton and hadron 
configurations of the colliding system which can contribute to the production of the inclusive particle 
with the momentum p. According to statistical physics, entropy of a system is given by a number 
of its statistical states as follows 



5' = In W^. (23) 

The most likely microscopic configuration of the system is given by the maximal value of 5. In our 
approach, the configurations comprise all mutually independent constituent subprocesses (2). The 
corresponding subprocesses which can contribute to the production of the inclusive particle with 
the 4- momentum p are considered as subject to the condition (3). The underlying subprocess, in 
terms of which the variable z is determined, is singled out from the corresponding subprocesses by 
the principle of maximal entropy S. 

In thermodynamics, entropy for an ideal gas is determined by the formula 

S = cvlnT + R\nV + So, (24) 

where cy is a specific heat and i? is a universal constant. The temperature T and the volume V 
characterize a state of the system. Using (22) and (23), wc can write the relation 

S = c ln[dN,h/dv\o]+ln[{l-x,)'^{l-X2y'{l-yar-{l-ybr']. (25) 

Next we exploit analogy between (24) and (25) . The analogy is supported by the plausible idea that 
interaction of the extended objects like hadrons and nuclei can be treated at high energies as a set of 
independent collisions of their constituents. Such concept justifies a division of the system into the 
part comprising the selected subprocess which imderlies production of the inclusive particle and the 
rest of the system containing all remaining microscopic configurations which lead to the produced 
multiplicity. The entropy (25) refers to the rest of the system. The multiplicity density dNch/dr]\o of 
particles in the central interaction region characterizes a "temperature" created in the system [36]. 
Provided that the system is in a local equilibrium, there exists a simple relation dNch/dri\a ^ T'' 
for high temperatures and small chemical potentials. Using the mentioned analogy, the parameter 
c plays a role of a "specific heat" of the produced matter. The second term in (25) depends on 
the volume of the rest of the system in the space of the momentum fractions {xi,x2, ya,yb}- The 
volume is a product of the complements of the fractions with the exponents which are generally 
fractional numbers, V = ll^ l^" ll'' ■ This analogy emphasizes once more the interpretation of the 
parameters ^i, S2, Ca, and Cb as fractal dimensions. In accordance with common arguments, the 
entropy (25) increases with the multiplicity density and decreases with the increasing resolution 0~^. 
In the considered analogy, the principle of a minimal resolution f2~^ with respect to all subprocesses 
satisfying the condition (3) is equivalent to the principle of a maximal entropy S of the rest of the 
colliding system. 

The entropy is determined up to an arbitrary constant Sq = In Wq ■ Dimensional units entering 
the definition of the entropy can be included within this constant. In particular, it allows us to 
account for a relation between the dimensionless multiplicity density dNch/dr]\o and the temperature 
T. This degree of freedom is related to the transformation 
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z^Woz, V ^ V'- 



(26) 



In such a way the scahng variable and the scahng function are determined up to an arbitrary 
multiplicative constant. The constant Wo is related to an absolute number of the microscopic 
states of the system. Its value is restricted by the positiveness of the entropy above some scale 
characterized by a maximal resolution fi"^. For the resolution corresponding to the fractal limit, 
Wo is infinity. Thus, the transformation (26) is connected with a renormalization of the fractal 
measure z in agreement with its physical interpretation (21). The renormalization represents a shift 
of the entropy by the corresponding constant, 

S = lnW + \nWo. (27) 



B. Chemical potential and phase space occupancy 

Let us consider a system where the number of particles can change. According to the first law 
of thermodynamics the internal energy U of such system can be expressed as follows 

dU = TdS - pdV + (TadNa + UadNa- (28) 

This shows how U depends on independent variations of the entropy S, volume V , and number of 
particles (antiparticles) Na, {Ng)- The parameters T and p are the temperature and pressure of the 
system. The potentials cja and cja consist of two parts [37] 

<Ja = lia+T\u{^a), (29) 

which are expressed via the chemical potentials jjba = —Us, and the phase space occupancy 7a = 75 

of the system constituents. First one controls the net number of particles (Na — Na) arising from 
the particle and antiparticle difference. Second term regulates the number of particle-antiparticle 
pairs {Na + Na). In statistical hadro-chemistry the factors are related to the relative and absolute 
chemical equilibrium, respectively [37] . The chemical potentials determine the statistical parameter 
Ta known as fugacity 

Ta = exp{aa/T) = 7^ exp(M„/T), (30) 

which governs the yield of a corresponding particle type. 

As shown in Sec. IV, the scaling functions of different hadrons coincide each other when applying 
the transformation (17) with the appropriate values of air. The coefficient ap ratio of the constants 
Wo (27) for single hadrons. In principle, these properties can be used for estimation of the ratios of 
the corresponding phase space occupancy parameters ja- In Boltzmann approximation, the hadron 
yields are proportional to the chemical fugacities. If one assumes that at sufficiently high energy 
the chemical potentials /Za become negligible on the right-hand side of (29), then it is possible to 
write 

71/72 = aFi/aF2- (31) 

The ratio describes the relative yields of hadron pairs with the flavors Fi and F2. Exploiting the 
connection between the temperature T and the multiplicity density dN/dr]\Q discussed after Eq. 
(25), we can rewrite the relation (29) into the form 

(TFi = + g{dNch/diq\ofl^ \n{aFjaF^). (32) 

Here 5 is a constant which does not depend on the particle type. Knowing the values of ai? for 
different hadrons (they arc quoted in Figs. 2 and 3 relative to tt^ mesons with a^r = 1), one can 
obtain a hierarchy of the potentials up for the individual hadron species. A systematic analysis of 
ap for mesons and baryons, its dependence on the mass, flavor content, and spin requires further, 
more detailed study. 
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C. ^;-Scaling and phase transitions 



The self-similarity of the particle formation at low transverse momenta is a new specific feature 
which can give information on thcrmodynamical characteristics of the colliding system and clarify its 
behavior. For sufficiently low px the scaling variable becomes resolution-independent function of the 
multiplicity density dNch/dr] and the parameter c. Both quantities characterize medium produced 
in the high energy collisions. The multiplicity density of particles in the central interaction region 
can be connected with a "temperature" of the colliding system. The parameter c is interpreted 
as a "specific heat" of the produced medium. Its value c = 0.25 was found to be the same for 
different types of the inclusive hadrons. The estimated error is at the level of 10%. We expect 
that this parameter should change in consequence of a phase transition in the produced medium. 
The increase of the multiplicity (or energy) density should affect the medium characteristics in 
that case. A possible violation of the ^-scaling in the low-^; region connected with a change of the 
"specific heat" c might be considered as a manifestation of such phase transition. Measurements 
of the cross sections for hadrons with heavy quarks at higher energies allow to reach even lower 
values of z where the saturation of t{j{z) could in principle change. We expect that such a change of 
the observed behavior of the scaling function could give indications on new physical effects in this 
region. 

There are quantities in thermodynamics which are sensitive to phase transitions. A discontinuity 
of specific heat, compressibility factor and coefficient of thermal expansion characterizes second order 
phase transition of the thermodynamic system. The second derivative of the chemical potential 

d^a/dT^\p = -cp/T (33) 

is connected with a discontinuity of the specific heat in such a case. The quantities c in (25) and ap 
in (32) signify a specific heat and chemical potential, respectively. Their potential dependence on 
multiplicity density can be studied from suitable experimental data. A discontinuity of c associated 
with a jump in the second derivative of ap, especially at high dNch/dr]\o, could be considered 
as a signature of a second order phase transition in the medium produced in high energy hadron 
collisions. 

It is usually assumed that bulk of the produced matter at low-pr (low-z) consists of multitude of 
the strongly interacting constituents. There is no direct information on the type of the constituents. 
Even at high collision energies y^, the single constituent processes at low transverse momenta are 
experimentally invisible. Only indirect signatures can reveal the nature of the constituents. We 
suggest to use the analogs of the thcrmodynamical quantities such as the specific heat and chemical 
potentials which can be useful in the description of the produced matter. These quantities can be 
sensitive to phase transitions in the multiple particle systems. 

According to the phase transition theory and theory of the critical phenomena, an abrupt change 
in the specific capacity, coefficient of thermal expansion and compressibility with a small change 
in temperature near to a critical point is typical for a second-order phase transition [38^0]. The 
second-order phase transitions have a discontinuity in the first derivative of the entropy 

cv = dS/dln{T)\v. (34) 

In our case, such signature would correspond to a discontinuity or sharp growth of the parameter 
c near some critical point Zc. One can speculate that the value of c should become dependent on 
the transverse momentum px and rise sharply enough with the multiplicity density at low pr- It 
is usually considered that the increase of the multiplicity density is connected with a growth of 
the energy density. Existence of the internal structure of the constituents would allow to store the 
energy into their internal degrees of freedom which is connected with a rise of the specific heat. We 
consider therefore that study of the inclusive cross sections of hadron production in the low-pr region 
at high multiplicity densities dNch/dr] and collision energies i/s is most preferable in searching for 
signatures of phase transitions in hadron matter. 

VII. CONCLUSIONS 

Here we summarize the main results obtained in the paper. New properties of the 2:-scaling - 
the flavor independence and the saturation of ip{z) at low z, were established. We have studied 
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the spectra of the inchisive particles produced in proton-proton and proton- antiproton coUisions 
in ^-presentation. The experimental data on inclusive cross sections of different hadrons obtained 
at ISR, RHIC, and Tevatron were used in the present analysis. The data cover a wide range of 
the collision energies = 19 — 1960 GeV, the transverse momenta pT = 0.1 — 10 GcV/c, and 
the production angles Ocma = 3" — 90". The variable z depends on the parameters S, ep, and c. 
The parameters S and characterize structure of the colliding (anti)protons and fragmentation 
process in the final state, respectively. Their values are fixed by the energy, angular, and multiplicity 
independence of 'tjj{z) in the high-p^ part of the spectra. The third parameter c is interpreted as 
"specific heat" of the produced medium. The z-scaling is consistent with c — 0.25 and S = 0.5 for 
all types of the analyzed inclusive hadrons. The value of increases with the mass of the produced 
hadron. Estimated errors of the determination of the parameters do not exceed 10%. It was found 
that the parameters arc independent of kincmatical variables. 

It was shown that renormalization of the scaling variable z allows to reduce the scaling function 
of different hadrons to a single curve. The scaling function '4i{z) was established to be flavor 
independent in a vide range of z. It means that the shape of 'ip{z) is the same for different hadrons 
in this region. This includes hadrons with light and heavy quarks produced in high energy pp- and 
pp-coUisions. A saturation regime of the function ip{z) was observed for z < 0.1. The approximate 
constancy of 'ip{z) was demonstrated up to z ^ 10"'^ for charmonia and bottomia. We conclude 
that a violation of the asymptotic behavior of the scaling function in the region of low z could give 
information on a phase transition of the produced matter. This should apply especially for the 
events with high multiplicity densities in both pp/pp- and Td^d-coUisions at high energies. 

A connection between the ingredients of the variable z allowing thermodynamical interpretation 
and the thermodynamical quantities (entropy, temperature, specific heat, chemical potential) char- 
acterizing multi-particle systems was discussed. The " specific heat" c was suggested as a parameter 
useful for characterization of these systems near a phase transition point. A microscopic scenario 
of the constituent interactions in terms of the momentum fractions was presented. 

The universality of description of the spectra of different hadrons over a wide kincmatical re- 
gion by the scaling function ip{z) reflects general features (symmetries) of the underlying physics 
phenomena over a wide scale range. Both soft and hard regimes manifest self-similarity of particle 
production at a level of hadron constituents. They are preferable to study collective processes and 
constituent substructure at low and high z, respectively. 

The obtained results may be exploited to search for and study of new physics phenomena in 
particle production in the high energy proton-proton and proton-antiproton collisions at U70, RHIC, 
Tevatron, and LHC. 
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